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Abstr act. We investigate in details the properties of the wavelengths (A = 10 — 20 /im) recently followed by Schnci 



disk su rrounding the 8Myr old star HR4796A, one of the 
few stars bringing precious clues to better understand the 
scenario which lead to planetary system formation. We 
propose a model able to reproduce all the available ob- 
servations : the full spectral energy distribution from the 
mid-infrared to the millimeter wavelengths, resolved scat- 
tered light and thermal emission observations. We show 
that the circumstellar matter splits into two dust compo- 
nents : a cold annulus, peaked at 70 AU from the star, 
made of ISM-like grains (amorphous composition, poro- 
sity ~ 0.6) larger than 10 /im and a population of hot dust 
close to the star (at about 9 AU) made of comet-like grains 
(crystalline composition, porosity ~ 0.97). Both dust pop- 
ulations are highly collisional and the grain size distribu- 
tion in the cold annulus is found to be cut-off by radiation 
pressure. At 70 AU, bodies as large as a few meters are re- 
quired to fit the data leading to a minimum disk mass of a 
few Earth masses and to a gas to dust ratio less than 1. We 
discuss aftewards some implications on the disk structure 
and effects of larger bodies. 



dcr ct al. (1999 ) in scattered light (A = 1.1 and 1.6 /im). 
Interestingly, the disk lies nearly in the direction to 
HR4796B, a phy sically bound companion of HR4796A 
( Jura ct al. 1993J ), located 7.7" away, i.e. 515 AU in pro- 
jected distance according to the Hipparcos star distance 
d* = 67.ljj;5pc. 

According to the resolved images, the disk is extending 
only a few tens of AU outside the peak of the dust dis- 
tribution assessed to about 70 AU ( Schneider et al. 1999| ). 
Outer truncation of circumstellar disks within a binary 
system has been predicted at typical distances of 1/3 - 



1/2 of the binary separation (Papaloizou & Pringle 1977; 



Artymowicz et al. 1994). However, the impact of the com- 



panion HR4796B on the disk extension is still unclear. 

In the inner part of the disk, Jura et al. (1995| ) first 
suggested a depletion of dust close to the star so as to re- 
produce the 110 K color temperature deduced from IRAS 
data. The comparison between resolved mid-IR images 
and first order modeling (see below) confirmed the neces- 
sity of an inner hole in the disk at 55±15 AU from the star 
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HR4796X 



for Koerner et al. (1998) and 60±20 AU for Jayawardhana 



According to Koerner et al. (1998), a second popula- 



1. Introduction 

HR 4796 A is one of the A type stars showing a high in- 
frared excess (£disk / L* = 5.10~ 3 ) due to the thermal 
emission from circumstellar dust. The excess is twice that 
of f3 Pictoris, the best studied ca se among Vega-like s tars 
so far. The star age, 8±2Myr flStauffcr ct al. 1995|) , is 
about 5 to 15 times less than (3 Pictoris implying that 
HR 4796 A is tracing an evolutionary phasis prior to that 
of (3 Pictoris. 



Koerner et al. (1998) and Jayawardhana et al. (1998 ) 
resolved first the circumstellar dust at thermal infrared 



tion of hotter grains may however lie closer to the star 
(inside the inner hole). Located at distances similar to 
those of the zodiacal dust in our Solar System, this dust 
population would be responsible for both the 12.5 /im de- 
tected excess and for a faint emission at 20.8 /im (less than 
10% of the total flux at the wavelenght) in excess of the 
resolved disk and centered on the star. 

The grains properties are so far poorly constrained. 
According to Jura ct al. (1995| ) , grains smaller than about 
3 lira, are blown outward by radiation pressure. Consid- 
ering the Poynting-Robertson effect, they found that the 
grains are probably larger than 40 /im under the assump- 
tion of a 40 AU inner hole. Koerner ct al. (1998| ) used 
thermal absorption/emission laws for the grains following 
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the model proposed by Backman et al. (1992) for /3 Picto- 
ris : Qabs is constant for A < Ao and proportionnal to 
(A/Ao) otherwise, the parameter Ao is expected to be 
related to an effectif grain size ao in the disk. Although 
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Koerner et al. (1998) constrain Aq to 25±15/im, they do and the IRTF bolometer datum obtained by Jura et al 



not sp jcify the ratio ^ in the case of HR4796A. This (1993 ) at 20 /im. They also include the available submil 



ratio could vary by a factor 40 between the lower and the limeter data : JCMT upper limit at 800 /im (Jura et al 



higher value (Backman et al. 1992|). 



Before beeing imaged in the mid-IR spectral range, no 
optical or ncar-IR observations succesfully resolved the 
disk around HR4796A, if we except for a slight asym- 
metry in the coronographic adaptive optics system im- 
ages per formed with the ESQ a daptive optics system in 
K' band Mouillct et al. (1997a ). Because of the low sig- 
nal to noise ratio in this data, these results were not pre- 
sented as a positive detection. Knowing the position angle 
of the disk, we reduced again these observations and we 
now marginally detect the disk. After a brief summary of 
thermal available data, we present in section the newly 
reduced scattered light images and compare then with 
more recent and better signal to noise ratio scattered ligth 



images (Schneider et al. 199G). 

As data become actually more numerous, we wish to 
better constrain the grain properties in the HR4796A 



disk. Sophisticat ed grains models developed by Green- 
berg et | al. (1972| ) succeeded both in repro ducing interstel- 



lar observations (extinction, polarization ( Li fc Greenberg 
1997|))]and in fitting the shape of the SEP of the disk 



surounding 3 Picto ris JLi fc Greenberg 1998 ; Pantin et al 



1997[jT|ln addition, Li fc Greenberg (1998J ) proposed a link 



between ISM particles and circumstellar evolved environ- 
ments assuming that dust gr ains in the 3 Pictoris disk ar e 
of cometary origin (see also Backman fc Paresce (1993|) ; 



Lecavelicr des Etangs et al. (1996] )) and could be fluffy 
aggregates of primitive interstellar dust. 

In this paper, we adopt a similar approach for the 
grains in the HR4796A disk. We describe section p| the 



4.1 



disk model assumptions and try to fit (sections 
4.2) the Spectral Energy Distribution (hereafter SED) to 
derive some physical and chemical grain properties such 
as typical size, porosity, presence of ice and finally to es- 
timate if the grains are more similar to interstellar dust 
grains or to comet-like grains. 

In addition, we test whether the presence of a second 
population responsible for the 10 /zm excesses is necessary 
or not. Given the constraints on the grain distribution 
derived from the SED fitting, we try to reproduce the 
thermal and scattered light resolved images (section |4.3| ). 
We finally discuss section || the implications on the disk 
dynamics. 

2. Available data 

2.1. Spectral Energy Distribution 

The available photometric data for the HR 4796 A disk are 
summarized in Table [fl. They include color-corrected IRAS 
measurements, IR TF bolometer measurements (|Fajardo- 
Acostact al. 1998J ), Ke ck/Mir lin (|Kocrner et al. 199S| ) and 
Cerro Tololo/ OSCIR ( Jayawardhana et al. 1998J ) images 



1995) and recent SCUBA measurements at 450 /xm and 
850 /im ( prcaves et al. 1999Q . 



A 

[/im] 



Excess 
*o bs [Jy] 



Uncertainty 
A$ obs [Jy] 



Reference 



4.80 


no excess 




7.80 


0.067 


0.037 


9.80 


no excess 




10.1 


0.087 


0.026 


10.3 


0.057 


0.024 


*11.6 


0.086 


0.070 


*12.0 


0.173 


0.028 


*12.5 


0.101 


0.018 


*12.5 


0.133 


0.027 


*18.2 


1.100 


0.150 


20.0 


1.860 


unknown 


*20.8 


1.813 


0.170 


*24.5 


2.237 


0.700 


*25.0 


3.250 


0.130 


*60.0 


8.630 


0.430 


100.0 


4.300 


0.340 


450.0 


0.180 


0.150 


800.0 


<0.028 




850.0 


0.0191 


0.0034 



Fajardo-A. 
Fajardo-A. 
Fajardo-A. 
Fajardo-A. 
Fajardo-A. 
Fajardo-A. 



et al. 

et al. 

et al. 

et al. 

et al. 

et al. 



(1998) 
(1998) 
(1998) 
(1998) 
(1998) 
(1998) 



IRAS 



Koerner et al. (1998 ) 
Fajardo-A. et al. (1998) 



Jaya wardhana et al. ( 1998 ) 



Jura et al. (19931) 



Koerner ct al. (1996) 



Koerner et al. (1996) 



IRAS 
IRAS 
IRAS 



Greaves et al. (19991) 



Jura et al. (1995Q 



Greaves et al. (1999) 



Table 1. Available infrared and submillimeter measure- 



ments. Refer to section 4.1 for the meaning of the mark 



anc ] 2.2. Resolved mid infrared thermal data 



Mid-IR images revealed the morphology of the dust dis- 
tribution sensitive to these wavelengths. Actually, the 
rather low resolution does not allow to precisely describe 
the radial shape of the surface brightness along the ma- 
jor axis and its vertical shape. However, simple models 
are consistent with data if the suface density of grains 
follows a radial power law r ~ r where T is badly con- 
straincd in the range [0, 2. 5[ ( Koerner ct al. (1998 ) and 
[layawardhana et al. (1998 )). Assuming that the disk is 
optically and geometricall y thin (2-D disk mode l) and us- 
ing a Bayesian approach, Koerner et al. (1998 ) assessed 
a disk inclination less than about 20° with respect to 
the line of sight. Indeed, the vertical height to radial 
height ratio for a given isophote at 20.8 /im is close to 
0.35 ± 0.05 which corresponds to a flat disk inclined at 
sin- 1 (0.35±0.05) ~ 20.5 ±3° from edge-on. Finally, both 
Koerner ct al. (1998 ) and layawardhana et al. (199q ) im- 



ages are consistent with a disk located at the position an- 
gle (PA) in the ra nges : [22°, 34°] for [Koerner et al. (1998| ) 
and [20°, 40°] for [layawardhana et al. (1998|). 
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2.3. Resolved near infrared data 

2.3.1. Reduction procedure 

The coronographic observations we reinvestigate in this 
paper have been first presented in Mouillet et al. (1997a|) . 



This paper also describes the reduction procedure used. A 
main step in coronographic data reduction is the removal 
of the scattered light remaining around the mask. This 
implies to use a comparison star which scattered light is 
scaled to that of the object of interest. The scaling factor is 
generally estimated by azimuthally averaging the division 
of the star of interest by the comparison star. Azimuthal 
averaging over 360° prevents from detecting pole on disks 
and decreases the chance of detection of inclined disks or 
faint structures. Knowing the orientation of the disk, we 
were able to estimate more precisely the scaling factor in 
a section perpendicular to the disk. 

2.3.2. Radial excess almost in the direction of HR4796B 

We present Figure |l| the newly reduced coronographic im- 
age of HR4796 A in K' band. An excess is present in the 
NE and SW regions. Due to the spreading of the ex- 
cess and due to the rather low signal to noise, the po- 
sition angle is badly constrained but is roughly consis- 
tent with both the PA deduced from previous images 
and the direction of the physical companion HR4796B. 
We extract the residual signal azimuthally averaged in 
the ranges of PA : [20°, 40°] and [200°, 220°] as a func- 




Fig. 1. Newly reduced coronographic image in logarithmic 
scale of HR4796A in K' filter revealing an excess close 
to the expected PA. The field-of-view is 13" x 13" with a 
sampling of 0.05" per pixel, smoothed to 0.2". The nu- 
merical mask occults up to about 60 AU (0.9") in radius. 
One can note the faint companion HR4796D (K'=14.5) at 
PA=311° detected by |Kalas k Jcwitt (1993[) and pVlouillct 



tion ofl the distance from the star, We plot in Figures B 
and @ the comparison with the noise on the image es- 
timated in the regions perpendicular to the excess. Be- 
tween about 1.1" and 2.3", the NE and SW remaining 
signals appear above the noise whereas in the perpendicu- 
lar region the residual signal fluctuates around which 
characterizes the lack of circumstellar structure in the 
NW and SE regions. Although the signal to noise ratio 
is rather low close to the star, we find an azimuthally 
averaged surface brightness at 1.1" - 1.2" (75 - 80 AU) 
from the star : 12.2±0.5mag.arcsec -2 quite similar to the 
range 11.7 - 12.7 mag . arcsec~ 2 measured in H band by 
Schneider et al. (1999| ) between 1.1" and 1.2". The ob- 



et al. (1997a) 



dimensional distributions and no the grain density at the 
normalized distance r = ro when z = 0. The distance ro 
is fixed to 70 AU in all the simulations and R(r) is nor- 
malized so as to have R(ro) = 1. 

Following Artymowicz et al. (1989 ), we assume a ver- 
tical distribution characterized by an exponential with a 
shape depending on the parameter 7 (7 > 0) : 



Z(r, z) = exp 



served surface brightness is roughly proportionnal to r 
but deconvolved i mages would lead to steeper slopes as 
also suggested by Bchncidcr et al. (1999 ) observations. 



C(r) 



with C(r) = Cro - 







3. Disk model assumptions 

3.1. Grain distribution 

We adopt a parametrical approach to describe the grain 
distribution in the disk similarly to that adopted to de- 
scribe the /3Pictoris disk ( Artymowicz et al. 1989J ) : 

n(r, z) = no R{r) Z{r, z) 

where n(r, z) is axisymmetrical grain density, R(r) and 
Z(r, z) are respectively the radial and the vertical no- 



where the radial shape of the vertical scale height £( r ) 
depends on the parameter (3 > 0. ( ro is the normalized 
height when r — r^. This notation is useful because the 
integration of Z(r, z) along the vertical axis (z) of the 
disk is proportionnal to r" . Therefore, the surface density 
of the grains, which is the interesting quantity to model 
the photometric data, can be written as : 



cr(r) 



+ OC 



n(r, z)dz — a$ x R(r) — 
Jo 



where ctq = C 7 noCr an d C 7 = 2 J exp(— x 1 )dx (x is 
an non-dimensional variable). 
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1.5 2.0 

Distance from the star [arcsec] 

Fig. 2. NE residual signal in the region around PA=30° 
compared to the noise (plain line). The latter is inter- 
preted as the detection limit of the disk. The residual sig- 
nal in the regions perpendicular to the disk is close to 
on average (not represented in this logarithmic scale). 
On contrary, the residual signal in the NE region is con- 
tinuously greater than the noise between 1.1"— 1.2" and 
2. 3" -2.4". 




1.0 



1.5 

Distance 



2.0 
from the star [arcsec] 



2.5 



Fig. 3. Same as Figure but for the SW region (around 
PA=210°). 



sublimation temperature (~ 1800 K). The maximum of 
the surface density cr(r) is not r c but also depends on a[ n , 
a out and /3 : 



^max(cr) 



-r„ 



(2r in -2r„ u t)- 



(1) 



with Tin = a m + (3 and r out = a ut + P- Outside this max- 
imum, a(r) is roughly proportionnal to r r ° ut , and roughly 
proportionnal to r Fin inside. Finally, the SED computa- 
tion requires a four parameters grain distribution : r c , Tin, 
r out , and er . 

3.2. Grains characteristics 

3.2.1. Composition assumptions : 

We assume that the grains are porous aggregates made 
of a silicate core coated by an organic refractory mantle 



(Greenberg (1986) and references therein, Pollack et al 



(1994)). Vaccum is assumed to fill the holes due to poros- 



ity. If the grain temperature is less than the sublimation 
temperature of the H2O ice (T su b,H 2 o — 110 — 120K), the 
H2O ice is able to fill all or part of the vaccum in the grain. 
The relative abundances of each compound are given by 
the volume ratios : gsior = ^ and gsiorH 2 = V Vh^' ' 
where Vgj, V OI and Vh 2 o are respectively the silicate, 
the organic refractories and the H2O volumes. The cor- 
responding grain density p g is : 



: Pg = (1 - p) 



ffSior PSi 
1 + gSior 



Pa 



1 + qsh 



Ph 2 o 

<?SiorH 2 



(2) 



p Si = 3.5g.cm '\ 
" 3 . Given the Si+or 



where the individual densities are : 
por = 1.8g.cm -3 and pn 2 o = 1.2g.cm" 
aggregate porosity P, the corresponding percentage of vac- 
cum removed due to the presence of the H2O ice is there- 
fore : 



PH 2 o = 



- 1 



100 



<7SiorH 2 



% 



.Li & Greenberg (iyys) proposed two extreme types of 
dust in Vega-like disks : amorphous or crystalline grains 
mainly depending on the evolution of the circumstellar 
environment and probably depending on the spectral 
type. In this paper, we investigate the two following types 
of grains : 



We parametrize the radial distribution of the grains 
with a smooth combination of two power laws : 



R(r) ex - 




with a in > and a out < 0. The inner disk is assumed 
to be cut-off at the distance corresponding to the grain 



1. "ISM-like grains" : 

amorphous Si and amorphous H2O ice, 
for spherical grains and qsior 



Qsic 



j_ 
2.12 



1 

2.')] 



for 



cylindrical grains (Li & Greenberg 1997), 

— porosity P between 0.4 5 ( fMathis 1996| ) and 0.8 
( fVIathis fc Whiffen 1989| ). 

"comet-like grains" : 

— crystalline Si and crysta lline H2O ice, 

— gsior — \ (comet Halley, JGrccnbcrg (1998 )), 
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Fig. 4. Grain temperature versus the distance from the 
star for amorphous and crystalline grains and for different 
grain sizes a. The grain temperature has been computed 
from the energetic balance between the absorbed flux and 
the emitted flux over all the spectrum. We use a Kurucz 
spectrum for an A1V star scaled to the V magnitude and 
extrapolated with a A -4 power law (Rayleigh- Jeans) to 
simulate the star flux. The break in the grain temperature 
is due to H2O ice sublimation (see also Figure pi). The dust 
temperature of the icy grains follows a radial power law r q 
with q between —0.36 and —0.47 for the amorphous grains 
and q between —0.40 and —0.43 for the crystalline grains. 



porosity P between 0.93 and 0.975 (comet Halley, 
Greenberg fc Hage (199q )). 



For all the models, we consider spherical grains and there- 
fore we fixed the volume ratio gsior to i. 



3.2.2. Optical properties computation : 

We compute theorical absorption/emission and scattering 
efficiencies for spherical grains via Mie theory (Bohren & 
Huffman 1983| ) knowing the ratio x = f where a is the 
grain radius and knowing m(A) the complex index of re- 
fraction of the grain. In some extreme cases, Mie theory 
is not efficient : we used the Rayleigh- Gans theory when 
\m(X)\x > 1000 and |m(A) - l\x < 0.001 or the geomet- 
ric optics when |m(A )|a: > 1000 and \m(X) - l\x > 0.001 
( [Laor fc Drainc 1993J ). 

We adopt the Maxwell-Garnett effective medium 
(MGEM) theory to compute m(A) for the mixture (Si + 
"or" + H 2 ice and/or vaccum) given the complex indexes 
of refraction of each compound, the volumes ratios and the 
porosity. We prefer the MGEM theory to the Bruggman 
effective medium (BEM) theory because the MGEM the- 
ory considers that a mixture is composed of an inclusion 
(the core) embedded in an homogeneous matrix (the man- 
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10% of ice 
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Fig. 5. Grain size corresponding to the H2O ice sublima- 
tion limit versus the distance from the star for a typical 
amorphous grain and a typical crystalline grain. On the 
left side of the curve, the ice is sublimated; on the right, 
grains are icy. 



tie) whereas the BEM theory treats both components in 
the same way (Bohren & Huffman 1983). 



3.2.3. Optical constants : 

The wavelength dependence of the grain complex index of 
refraction m(A) depends on the individual compounds : 

— > Silicates : there is presently no consensus on what 
could be "representative" astronomical silicates. In the 
case of amorphous silicate s, we arbitrarily adopt the modi- 
ficd (Laor fc Drainc 1993]) Draine fc Lee optical constants 
( Drainc fc Lee 1984 ; Draine 1985 ). In the case of crys- 



talline silicates, we adopt the analytic complex indexes 
of refraction for crystalline olivine particles proposed by 
Mukai fc Koike (1990| ) be tween A = 7 ^m and A = 200 ^m 
and we extrapolate (as Li fc Greenberg (199§| )) with 



Drainc fc Lee (1984 ) optical constants for A < 7/im and 
A > 200 /im. 

— ► Organic refractories : we adopt the optical constants 
computed by Li & Greenberg (1997) (a full discussion on 
organic refractories can be found in this paper). 

— ► Hi O ice : in the case of amorphous ice and between 
A = 1.3 /im and A = 144 ^m, we use the optical constants 
derived from |Fink fc Sill (1982fl (T = 82 K ) combined with 
the 100 K data from Hudgins ct al. (1993] ). In the case of 
crystalline ice and between A = 0.89 (im and A = 333 /jm, 
we adopt the hexagonal water ice optical constants at a 



temperature of about 100K derived by Ockman (1958) 



and Bertie et al. (1969). For both th e amorphous a nd crys- 



talline ices, we extrapolate with the Warren (1984 ) optical 
constants. 
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3.2.4. Grain size distribution : 

The choice of the grain size distribution n(a) in the disk 
is a key step in the grain modeling. We assume either 
a single grain size (a = ao) or a grain size distribution 
following the classical power law n(a)da ex a~ K da in the 
range [a m tn> a m ax] which may be better suited if the grains 
result from collisions among larger bodies. We fixed k = 
3.5, a typical value usually assumed for collisional disks 
flHellyer 1970[ pVIathis et al. 1977[ ). We use : n(a)da = 
2.5x _35 da; with x = a/a m i n so as to have : 



Our aim is to minimize the \ 2 



2.5x _3 ' 5 da; = 1 with 



> 1 



Finally, the grains are characterized by their amorphous 
or crystalline nature and by the three free parameters : 
ao, P and ph 2 o for a single grain size distribution or by 
the four free parameters : a m - m , a max , P and pn 2 o for a 
—3.5 power law grain size distribution. Some examples of 
the dust temperature and of the grain size corresponding 
to the H2O ice sublimation limit versus the distance from 
the star are plotted Figures [| and || respectively. 

4. Modelling 

4-1- SED fitting assuming a single grain size 
distribution : 



Our aim is to derive the dust distribution and to constrain 
the grains properties with the available measurements. 

We already have a constraint on the peak of the grain 
distribution (r c ). Indeed, the shape of the isop hotcs in 
scattered light at A = 1.1 /im and A = 1.6 /imm ( Schnei- 
der et al. 1999| ) suggests faint anisotropic scattering prop- 
erties for the grains at these wavelengths (see also 4.3.2 ) 
implying that the observed maximum surface brightness 
at roughly 70 AU is probably close to the peak of the grain 
distribution. We therefore fix r c = 70 AU. 

We let the 6 remaining parameters free in their respec- 
tive ranges : 



— the 2 grains distribution parameters : 

r in and r out in [2 , 13] and [—13 , —2] respectively, 

— the 3 grains characteristics parameters : 

grain size ao in [1 /im , 600 /im] , porosity P in the range 
[0.45 , 0.8] if the grains are amorphous or in [0.93 , 
0.975] if the grains are crystalline and percentage of 
H 2 icepH 2 o in [0% , 50%]. 

— the normalized dust surface density ao : 

given a set of the 5 previous parameters, the normal- 
ized dust surface density ao is used as a flux scaling 
factor so as to provide the best fit of the SED. 



As a first step, we try to fit the 10 - 100 /im SED only. 



X 



N 

E 

»=1 



$obs(A t ) - $ sim (A I 
A$ obs (A;) 



where $ s im(Ai) are the simulated flux and N the number of 
data and discuss afterwards the consistency with thermal 
and scattered light resolved images. 

We start the minimization with a random vector of 
the parameters. One has to note that we limit A?" to the 
ten data marked with a "*" . Indeed, taking into account 
the numerous measurements at about 10 /im, compared to 
the few data at larger wavelengths, would give a too high 
weight to this SED region. The x 2 probability p(x 2 , xV— 6) 
gives us a criterion for the quality of the fit. For example, 
so as to have p > 10% we need \ 2 < 7.78 with N = 10. 
We present Table |] six different sets of the parameters 



0.00 


— * 

= + 


d 



:ta used to derive the parameters 

her data ,--<i"~ - " """""•^ 


**- 


1.00 


J* I 

- f ;i 0.10 

- / 1 

: / 

_ / 0.01 




- 


0.10 

01 


r I 

j /■, i 


-py 9 10 11 12 13 u 

/ /* Wavelength [/J.m] — 

model #1 (amorphous grains) - 

model §4 (crystalline grains) 



10 



Wavelength [/^m] 



100 



Fig. 6. Two fits of the 10 /im - 100 /xm Spectral Energy 
Distribution of the HR 4796 A disk assuming a single grain 
size distribution a — ao and a single grain distribution 
peaked at r c = 70 AU from the star (parameters from 
Table |). 



which lead to the best 10 - 100 /xm SED fitting assum- 
ing a single grain size distribution ao and a single grain 
distribution. 

The 3 crystalline grains models satisfy the \ 2 crite- 
rion : p > 10 % whereas the amorphous grains models do 
not. Amorphous grains do not provide a proper fit of the 
10 - 12.5 /xm data. The crystalline hypothesis reproduce 
very well the 10 /xm emission feature. Furthermore, the 
10.1 /xm and 10.3 /xm data, not used to derive the param- 
eters, are also well fitted with the model (Figure ^|). Note 
that the 20 /im - 100 /xm continuum is fitted both with 
crystalline grains and with amorphous grains. 

Grains with Oq S > 60 /xm are needed to reproduce the 
shape of the SED (see Table g for the definition of ag ff ) 
which leads to grains sizes ao larger than 150 /xm. Most 
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#5 
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#6 
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15 


7.5 
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0.23 


4.09 



Table 2. Sets of parameters providing the best fits of the 10 - 100 pm SED assuming a single grain size distribution 
and amorphous grains (models #1, #2 and #3) or crystalline grains (models #4, #5 and #6). We also add some 
interesting quantitites such as : the effective grain size a eff = a (1 — P) 1 ' 3 which represents the size of a not porous grain 
with the same mass as the grain with the size ao and the porosity P, the final grain porosity P W ithH 2 o = P (1— Ph 2 o/100) 
taking into account that ph 2 o % of the vaccum has been removed and filled with H2O ice, the grain density p g (equation 



||) and the disk mass M<j(ao). For all the models, the peak of the surface density is close to 



K(a) 



= 67AU. 
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Table 3. Sets of parameters providing the best fits of the full SED assuming a collisionnal grain size distribution 
(oc a~ 3 - 5 ) and amorphous grains (models #7, #8 and #9) or crystalline grains (models #10, #11 and #12). Refer 
to Table gfor the definition of the quantities a off , P W ithH 2 o and p g . Md(amin) is the mass corresponding to the smallest 
grains and Md = 2Md (a m in) \7a m ax/amin is the total disk mass. The peak of the surface density r mm („) is close to r c . 



of the grains are located further than the H2O ice sub- 
limation distance (Figure 0) and have to be significantly 
icy : j»h 2 o > 15 — 20% (consistent with Jura et al. (1998) 
assumptions). 

We find Tin to be close to 3 and the outer surface 
density (r out ) to be steeper than r~ 7 (more probably as 
steep as r -10 ). Whereas the outer slope of the grain sur- 
face d qnoity io conoiotont with the ocattorod ligth imagoa 



the in ner elopo is not. It would induce much more light : 3' 



the star. Under these assumptions neither the crystalline 
nor the amorphous grains provide good fits of the 10 - 
100 /mi SED (Figured). In the present case, we also ex- 
clude to add an inner dust distribution (typically inside 
50 AU) to fit the SED at short wavelengths. In this scheme, 
the inner and outer disks would have similar contributions 
(flux) at 20 /im (Figur e jS|), which is not compatible with 
20 jum reso lved images ( Kocrncr et al. 1998 ; Jayawardhana 



in the inner disk than observed. Indeed, an r 3 inner sur- 
face density yields an inner midplane surface brightness 
in scattered light less steep than at least r 2 ~P (case of an 
edge-on dis k assuming iso tropic scattering properties for 
the grains, Nakano (1990J )). As an example, we plot Fig- 
ure [7] the surface brigthness contours for a 15° tilted (from 
edge-on) disk in scattered light at A = 1.1 pm. 

Taking into account the constraint from scattered light 
images, we therefore try to derive the grains proper- 
ties which could reproduce the 10 - 100 pm SED with 
Tin ~ 10. Such an inner slope acts as a cut-off for the short 
wavelengths efficient emitters which arc located close to 



In conclusion, assuming a single grain size distribution, 
we are able to properly fit the 10 - 100 /im SED but not the 
resolved images. Moreover, the presence of hot crystalline 
grains seems necessary to fit the 10 pva measurements. 



4-2. SED fitting assuming a grain size distribution 
proportionnal to a -3 ' 5 : 

We now assume a collisional grain size distribution n(a) 
proportionnal to a -35 between a minimum and a maxi- 



mum gram size a n 



and 



We also assume that the 
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o.o 
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Fig. 7. Surface brightness contours of a 15° tilted disk 
in scattered light (A = 1.1 /xm) assuming the set of pa- 
rameters #4 (Table ||) and isotropic scattering properties 
for the grains. We also assume a flat disk : j3 = 0.75, 
£ ro = 4AU and an exponential vertical distribution : 
7 = 1. The contour levels are in an arbitrary unit ; the 
upper level is fixed to 100. In the inner part of the disk, 
the surface brightness decreases slowly with the distance 
from the star : for instance, the surface brightness at 0.5" 
from the distance corresponding to the upper level still 
represents 60% of the maximum surface brightness. Such 
an inner dust distribution (Tj n ~ 3) is then not consistent 
with HST observations ([Schneider et al. 1999). 



850 /xm measurement. We keep the porosity P and the 
percentage of H2O ice free in the same ranges as for the 
single grain size distribution and investigate the respective 
broad ranges [1 /im , 500 /xm] and [600 /im , 5 m] for the 
minimum a m i n and the maximum a max grain sizes. 

The results are summarized in Table [3| Whatever the 
nature of the grains, the fits do not satisfy the % 2 crite- 
rion (p > 10% if x 2 < 9.24), but micronic to centimet- 
ric amorphous grains better reproduce the shape of the 
SED continuum than crystalline grains (Figure ||) . There- 
fore, crystalline grains (if any) probably represent a mi- 
nor part of the grain population in the annulus peaked at 
70 AU. On the other hand, the 10 /xm measurements are 
not well reproduced with the amorphous grain model still 
suggesting (see previous section) the presence of a small 
amount of crystalline grains in order to significantly emit 
at 10 /xm but not at larger wavelengths. If there are crys- 
talline grains at 70 AU, they are neither numerous nor 
hot enough to produce such a spectral energy distribu- 
tion. Therefore, if the 10 /im data are not overestimated 
(especially the 10.1 /xm, the 10.3/xm and the 12 /xm mea- 
surements), we can exclude to fit the full SED from the 
NIR to the submillimeter wavelengths with the same dust 
population. 
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Fig. 8. Same as Figure g but assuming T m = 10. 



relative abundances of the different compounds and the 
porosity do not depend on the grain size. 

We adopt a spatial grain distribution consistent with 
the scattered light images and therefore fix r c to 70 AU, 
the inner slope T m to 10 and the outer slope r out to —11.5. 
The grain properties are now the free parameters left. 

4.2.1. Single grain population : 

We try to fit the SED by minimizing the \ 2 associated 
to the 10 data marked with a "*" in Table [y plus the 
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Fig. 9. Two fits of the full SED assuming a collisional 
grain size distribution and a single dust population peaked 
at 70 AU from the star. The 10 /xm and 850 /xm data are 
not properly fitted. 



4.2.2. Two grain populations : 

— > The cold annulus : 

We now assume that the excesses at large wavelengths 
mainly come from the annulus resolved in scattered light. 
This is consistent in terms of spatial extension with the 
disk observed at 20.8 /xm (Kocrner et al. 1998) suggesting 
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Table 4. Sets of parameters providing the best fits of the 20 - 850 /im SED assuming a collisionnal grain size 
distribution (ex a -3,5 ) and amorphous grains (models #13, #14 and #15) or crystalline grains (models #16, #17 
and #18). Refer to Table || for the definition of the quantities a cS , -P w ithH 2 0> Pg, Mi(amin) and Md. The peak of the 
surface density r max t a -\ is close to r c . 



that the same dust population may be responsible for both 
the 20.8/im and the 1.1 /im images. 

We therefore try to fit the 5 data at wavelengths larger 
than 20 /im (excluding the 24.5 /im and 450 /im data due 
to too high uncertainties) and try to constrain the grains 
properties consistent with the colder part of the SED. We 
fix the porosity P to 0.6 for amorphous grains or to 0.93 
for crystalline grains as suggested by the full SED fitting 
(Table 0). We keep the three parameters a m - m , a max and 
PH-.Q f ree in the same ranges as in the previous subsection 

The results of the x minimization are summarized in 
Table and Figure HG shows the best fits for each grain 
nature. The cold SED continuum is remarkably well fitted 
(p > 10% if x 2 < 2.71) with amorphous grains whereas 
the crystalline grains lead to very poor % 2 . Very little icy 
grains (ph 2 o < 5%) with sizes larger than 10 /im (15 /im 
if we assume P — 0.7) are required to fit the 20 - 25 /im 
data and bodies as large as a few meters are necessary to 
correctly reproduce the submillimcter excess. This leads 
to an outer disk mass larger than a few Earth masses. 



— > The hot dust population : 

At 20.8/im, the simulated cold annulus (model #13) con- 
tributes to about 90% of the measured excess at this wave- 
length. Therefore and conversely to the single grain size 
model, the addition of a second grain population would 
not modify the consistency between the simulated ther- 
mal images at 20.8 /im and the observations. 

We therefore subtract the simulated flux coming from 
the outer annulus (model #13) to the observed excesses 
and try to fit the remaining excesses with a second dust 
distribution made of crystalline grains (suggested by the 
single grain size distribution fitting (section 4.1)). The full 



SED assuming two grain populations is now very well fit- 
ted (see Figure O) . In this example, the inner dust distri- 
bution represents a mass of about 3.5 x 10 -5 M ffi made of 
450 /im crystalline grains (P = 0.97) peaked at 9 AU from 
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Fig. 10. Two fits of the 20 - 850 /im SED assuming the 
dust resolved in scattered light and peaked at 70 AU from 
the star is responsible for the excesses at wavelengths 
larger than about 20 /im. 



the star and heated at about 210 - 220 K (see for example 
Figure ||). 



Jf.,3. Images 

We now adopt t he tw o grain populations model from 
previous section (|4.2.2 ) : an amorphous dust population 
peaked at 70 AU (model #13) plus a crystalline dust pop- 
ulation at 9 AU from the star. We assume a disk inclina- 
tion itilt °f 15° with respect to the line of sight. 



4.3.1. Thermal images 

We simulate the disk in thermal emission so as to derive 
further informations on the vertical shape of the disk (pa- 
rameters ( ro , 7 and (3 assumed to be constant over all the 
disk). We actually have very few constraints on the ver- 
tical shape of the disk, also we arbitrarily start with pa- 
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Fig. 11. Left : the full SED fitting assuming two dust populations : model #13 + an inner dust population (see text). 
Right : simulated disks in thermal emission at A = 20.8 /im assuming grain properties and surface density derived 
from the SED fitting. We also assume a vertical structure quite simila r to the /3Pictoris dis k one. We simulate the 
disks with the same pixel size (0.14") as the observations performed by Koerner ct al. (1998 ). The images have been 
convolved with an Airy pattern to take into account the PSF of a 10 m telescope at A = 20.8 /xm. 
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Fig. 12. Simulations of the cold annulus peaked at 70 AU in scattered light at A = 1.1 /im for different asymmetry 
factor \g\ ; the inner hot dust not observable has not been ad ded. We assume the s ame spatial distribution and grains 
properties parameters as Figure [ll]. The pixel size is 0.076" ( [Schneider et al. 199S| ). 



rameters derived for j3 Pictoris disk : opening angle of the 
disk ~ 6° {i.e. C, ro ~ 4AU), (3 ~ 0.75 and 7 - 1 flMouil 
1997a| )- We plot (Figure [ll]) the simulated disk 



;1. 



let et 

at A = 20.8 /im. The simulated vertical height to radial 
height ratios for a given isophote are between 0.32 and 
0.35 ; they compare very well with the values 0.35±0.5 



from Koerner ct al. (1998) 



If we afterwards simulate the thermal images for dif- 
ferent sets of the parameters ( ro , f3 and 7, we get a few 
interesting results : 

— The normalized height ( ro at ro = 70 AU is proba- 
bly less than about 7 - 8AU leading to opening an- 
gles of the disk smaller than about 15°. For example, 
( ra = 10 AU provides a too large vertical structure and 
the vertical spreading of the dust yields a smaller sur- 



face flux density for the outer annulus compared to the 
inner one. Such a high constrast is not observed. 

— The parameter (3 is badly constrained : 0.75 provides a 
good shape of the isophotes (Figure |Tl|) but 1.5 might 
also be a correct value. 

— Given the resolution of the observations, the parameter 
7 is completly undetermined. 

4.3.2. Scattered light images 

To simulate the scattered light images, we adopt the 



empirical Henyey-Greenstein phase function (Henyey & 
jGreenstcin 1941 ) which depends on the single parameter 
\g\ (0 for isotropic scattering and 1 for total forward or 
backward scattering). Indeed, theoretical phase functions 
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from Mie theory are probably not realistic for porous fluffy 
grains. 

We simulate at A = 1.1 /im the cold annulus, which is 
responsible for the scattered light images, assuming dif- 
ferent values of the parameter \g\. As anticipated, due 
to the disk inclination with respect to the line of sight, 
the morphology of the observed disk strongly depends on 
the anisotropic scattering properti es (Figure fl2|) . The ob- 
served isophotes at A = 1.1 /mi (Schneider ct al. 1999) 
implies that the asymmetry factor \g\ is smaller than 0.15. 

At 1.1 /xm, the simulated cold annulus scattered a total 
flux close to 8.6 mJy (5.0 mJy at 1.6 /im) assuming \g\ ~ 0. 
The integrated flux density corresponding to the detected 
part of the annulus with HST (outside 0.65" in radius) 
represents 5.2 mJy at A = 1.1 /im, i.e. only slightly smaller 
than 7.5±0.5mJy from Schneider et al. (1999] ). 

We also simulate the disk in K' band. At 1.1" from the 
star along the major axis of the annulus, the flux density is 
about 16/iJy in a 0.05" x 0.05" pixel. The corresponding 
surface brightness is about 12.5 mag. arcsec -2 consistent 
with the 12. 2±0. 5 mag. arcsec -2 measured between 1.1" 
and 1.2" (subsection ^.3. 2j ). 



5.2. Radiation pressure 



4-4- Summary of model results 

Our model involves two dust components : 

— a cold annulus peaked at 70 AU from the star made 
of ISM-like grains in terms of chemical composition 
(amorphous) and porosity (P ~ 0.6) and a very few 
icy. Assuming a collisional grain size distribution (ex 
a -3 ' 5 ), we infer a minimum grain size in the ring of 
10 /mi and we find that bodies as large as at least a 
few meters may have already been formed. This leads 
to a total dust mass of a few Earth masses held in a 
very narrow annulus. 

— an inner dust population at about 9-10 AU from the 
star made of very porous (P ~ 0.97) crystalline grains 
responsible for the emission feature at A ~ 10 /im and 
for the s light excess centered on the star in the 20.8 /im 
images ( Koerner et al. 1998 ). 



5. Implications of the results on the disk 
dynamics 

5.1. Gas to Dust ratio 

It is of interest to assess the gas to dust ratio in the disk. 
Given the upper limit on the g as mass M g < 1 — 7M ffl 
derived by Greaves et al. (1999|) and the dust mass Mj ~ 
4M® from model #13 which c an be reasonably consid- 
ered as a lower limit (see also |5.4.1 ), we infer a gas to 
dust ratio less than 1. In this scheme, the gas would prob- 
ably not play anylonger an important role, as claimed by 
Greaves et al. (1999Q . 
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Fig. 13. Radiation pressure force to the gravitationnal 
force ratio (/3 pr ) assuming L* = 20 L© and M* = 2.5 M©. 
The amorphous grains characterize the cold dust annulus 
(solid ligne) and the crystalline ones the hot dust popula- 
tion (dashed ligne). Below the dotted ligne, the grains are 
gravitationally bound. 



The ratio /3 pr of radiation pressure to gravity as a func- 
tion of the grain size is given by : 



/3pr = 



3L*(Q P r(a)} 



with (Qpr(o)) = 



f Q P r(a)F4\)d\ 
fF*{X)d\ 



16nGM*cap g 

assuming L» = 20 L© and M* = 2.5 M©. We find a blow- 
out grain size limit (/3 pr = 1) close to 5 /mi for the amor- 
phous grains in the cold annulus (Figure |l3|). Also, amor- 
phous grains with sizes between 10 /mi (/3 pr ~ 0.5) and 
about 45 /im (/3 pr ~ 0.1) have orbits significantly excen- 
tric. Interestingly, the smallest grains we find in our fits in 
the cold annulus : a m i n ~ 10 /im is very close to the blow- 
out size limit. For the crystalline grains at 9AU from the 
star, the blow-out grain size limit is close to 65 /im (Fig- 
ure [TJh . The grains deduced from our modelling are then 
gravitationally bound. 

5.3. Times scales and dominant processes 

The life times of the grains in the disk under Poynting- 
Robertson drag, collisions and sublimation, can be written 
respectively as (e.g. Backman & Paresce 1993) : 



ipR 



4-7TC 



35- 



ih 



Vm r AU y r 



3 (Qfabs/i^* \Qabs, 

where (Q a bs) is the absorption efficiency averaged over the 
star spectrum, 

1 



ic 



>ii 



2.5 x 10 6 rl v 
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where f2 = yJGM*/r 3 / 2 is the Keplerian circular rotation 
frequency at radius r, and (a 2 ) the grain size square aver- 
aged over the grain size distribution. 
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where a eq = a (PpH 2 o/100) is the equivalent grain ra- 
dius for a pure ice sphere with the same mass of ice as in 
a grain with a size a. The time-scales for the cold annulus 
and the hot dust population are summarized in Table pi 
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Table 5. Time-scales under Poynting-Robertson drag, 
collisions and sublimation. At 70 AU : (a 2 ) ~ 5 a 2 ain , a eq — 
0.26 x a (model #13) and T g ~ 90 K (see Figure |). At 



9AU : Co — 0.5 grain. cm 



Og with oq — 450 /im, 



p g = 0.071 g.cm and H2O ice is sublimated (a cq = 0, see 
Figure ||). For both amorphous grains with a m i n = 10 /im 
and crystalline grains with a ~ 450 /im, we find (Q a bs) 
very close to 1. 



Collisions are then the dominant process in the cold an- 
nulus. The smallest grains produced by collisions in very 
short time-scales are not efficiently removed from the ring 
by Poynting-Robertson drag. Only radiation pressure is 
then able to explain the minimum grain size (10 /im) de- 
rived from the SED fitting. Also, this is consistent with 
the fact that this minimum grain size is very close to the 
blow-out size limit. Of course, this argument applies only 
if the gas does not play an important role, as suggested 
in a previous subsection (|5.1|). At 70 AU from the star, 
the H2O ice within the smallest grains is sublimated in a 
short time-scale in comparison with the star age. This is 
consistent with the little percentage o f H2O ice required 
to fit the 20 - 25 /im SED (subsection |4.2.2|) . 

The hot dust population is poorly constrained; the 
time-scales in Table ra have then to be taken with care. The 
life time of the 450 /im crystalline grains under Poynting- 
Robertson drag is long compared to the collision time- 
scale, and even 120 /im crystalline grains with /3 pr ~ 0.5 
have ipR about 4 times larger than t co \\. Also, collisions 
and radiation pressure are probably the two main pro- 
cesses which shape the hot dust distribution. 

5.4- Ring Dynamics and confinement 

5.4.1. Collision dynamics and Mass loss 

We here consider the implications of the model parame- 
ters for the ring centred at 70 AU given in Table 0. As 



representative example we consider model #13. This has 
a total grain surface density o- n = 1.79 x 10 3 grains. cm~ 2 . 
For the assumed size distribution (oc a -35 ), the mean col- 
lision time for the smallest grains with a = a m ; n — 10 /im 
is approximately t co \\(a m i u ) ~ 2.1 x 10 3 yr. As these grains 
are at the radiation pressure blow out limit and the ma- 
jority of grains have the minimum size, the minimum mass 
loss rate is given by : 



M 



M d (a m in) _ 47ra 3 Tlin p g n(a min ) 



t-coll (.^min J 



" ^coll v&minj 



(JqA 



where A = 1.7 x 10 30 cm 2 is the area of the annulus. For 
model #13 this corresponds to M ~ 2.5 x 10~ 6 Af e yr _1 . 
Thus the dust mass reservoir required to supply this mass 
loss over the age of the system is ~ 20M ffi , which corre- 
sponds to a surface density in solids E gr ~ 0.07 g.cm -2 . 

It is of interest to note that if we assume a gas to dust 
ratio of 100 in the primordial nebulae, this corresponds to 
a surface density E ~ 7 for the original protostellar disk. 
This surface density is consistent with that expected at 
70 AU for standard disk models corresponding to a mini- 
mum mass solar nebula (e.g. Papaloizou & Terquem 1999). 

If the surface density E gr resides predominantly in ob- 
jects with a = a max ~ lm, their collision time would be 
tcoll ~ 2ap g /(3f2E gr ) ~ 8.6 x 10 4 yr. This is still shorter 
than the age of the system. With a ring aspect ratio of 0.1 
for instance, collisions are expected to occur with a rela- 
tive velocity ~ 0.5 krn.s -1 , and are likely to be destructive. 
In order to have a collisional lifetime greater than the star 
age, the objects would have to have a larger than about 
100 m. A population of such larger objects providing a 
mass source for those with a~ lm, may exist but they 
would not be detectable. In that case the disk mass could 
be significantly larger than 4Af ffi . 

5.4.2. Formation, vertical height and possible presence of 
larger bodies 



In the standard model of planetary formation (Wcidcn- 



;chilling fc Cruzzi 1993| ), objects in the size range lm 



1 km are formed by accumulation of smaller objects 
through collisions when there was still gas present. The 
presence of gas drag, effective enough here with E ~ 
7 g.cm -2 , would be essential to damp eccentricity and in- 
clination so collisions occur with small relative velocity. 
However, the collision timescales would be the same as es- 
timated above. Thus it is reasonable that objects as large 
as 100 m could be formed if the gas residence time was 
similar to the present age of the system. 

When significant amount of gas was present during 
earlier phases, the disk would be expected to be thin be- 



cause of gas drag ( Weidenschilling & Cruzzi 1993) but 



this could have been thickenned due to gravitational per- 
turbations from large external objects. Unfortunately the 
observations are not yet able to provide contraints on the 
vertical thickness, but this does not affect the discussion 
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about mass loss rates or collision times because these are 
independent of the vertical thickness. 

5.4.3. Radial confinement 

The fact that the mass distribution in the ring is tightly 
radially confined to a width of about 0.2 times the mean 
radius, which may be comparable to its vertical thick- 
ness, gravitational confinement mechanisms may operate 
to sculpt its form. As it is likely that many collisions are 
required to reduce a metre sized object to fragments small 
enough to be blown out by radiation pressure, if the the 
vertical thickness is as much as 0.2r, the ring should be sig- 
nificantly wider without a confinement mechanism. Even 
if the ring is very thin as it might have been immediately 
post formation so that little collisional spreading occurs, a 
mechanism is still required to explain why the protostellar 
gaseous disk was truncated at around 70 AU. 

Confinement may occur through gravitational interac- 
tion with external perturbers. This situation is familiar in 



the case of planetary rings (Goldreich & Tremaine 1982) 



where the interaction with external orbiting satellites pro- 
vides confinement. It also occurs for accretion disks in 
close binary systems where interaction with a binary com- 
panion truncates the disk (Papaloizou & Pringlc 1977). 
The presence of a binary companion in the system dis- 
cussed here is suggestive that it could play a role in deter- 
mining the location of the outer boundary to the disk as 
well as pumping up inclinations after gas dispersal. Note 
also that an exterior as yet undetected giant planet might 
provide a similar role. 

The binary companion HR4796B is seen at a pro- 
jected distance of 515 AU from the central star whereas 
the outer disk radius is 70 AU. Thus if it is to play any role 
in confinement the orbit should be eccentric. As significant 
interaction with the ring requires the pericentre distance 
to be less than three times the ring radius, the minimum 
value for the eccentricity is thus e = 0.4, although it could 
be significantly larger on account of projection effects. 

Disk edge truncation can be effective in the presence of 
binary companions with significant mass ratio in circular 
or eccentric orbits. In the case of very eccentric orbits, as 
may well be the case here, the truncation process is most 
easily modelled as a succession of encounters which oc- 
cur at each pericentre pa ssage (pstrickcr 1994 ; Korycan 
sky & Papaloizou 1995; Hall 1997). From these studies 
it is clear that in order to be truncated, the disk radius 
should initially exceed about one third of pericentre dis- 
tance. Material beyond this radius is strongly perturbed 
and removed from the disk which develops a truncated 
edge there. The latter phenomenon, which occurs because 
of angular momentum transported from the disk to the 
companion, was seen in the gas disk simulations of Ko 



rycansky fc Papaloizou (1995 ) and the free particle simu- 
lations of Hall (1997 ) and is also expected for the disks of 
the type considered here. 



6. Conclusion 

We propose a two dust populations model able to re- 
produce all the available observations of the HR4796A 
circumstellar disk. The first dust population, peaked at 
70 AU from the star, shapes like a narrow annulus. Its 
small radial extension suggests a gravitational confine- 
ment mechanism both outside and inside the peak at 
70 AU. This ring is made of amorphous grains quite sim- 
ilar to those found in the interstellar medium. Under the 
assumption of a collisionnal grain size distribution, grains 
larger than 10 /um to meter sized bodies are required to 
fit the data. Collisions, occuring in short time scales, are 
expected to produce smallest grains down to the blow out 
size limit (~ 10/im). This conclusion applies because of 
the gas to dust ratio smaller than 1 inferred. The sec- 
ond dust population at 9-10 AU from the star is made of 
comet-like grains. Higher signal to noise ratio measure- 
ments and resolved images at A ~ 10 /im would allow to 
better constrain this dust population (spatial extension, 
grain size range) and its dynamics. 

The dust surrounding HR 4796 A probably shows an 
intermediate stage for the circumstellar matter between 
pre-main sequence stars and evolved ones like /3Pictoris, 
Fomalhaut, Vega, 55 Cnc . . . The presence of a companion 
(HR4796B) is also an interesting peculiarity among the 
Vega-like stars. As a large number of field stars are mem- 
bers of binary (or more) systems, HR4796 A does not rep- 
resent a marginal case. Therefore, the dynamics of the disk 
has to be further studied, in terms of chemical evolution of 
the grains and dust distribution to better understand the 
steps which lead or not to planetary systems formation. 
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